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CHAPTER  I 

INTRODUCTION 
I . Purpose 

Perhaps  the  most  common  configuration  of  vehicles 
capable  of  supersonic  flight  is  that  of  a supersonic  nozzle 
embedded  in  a blunt-based  afterbody,  as  typified  in  Fig.  1. 
The  complex  flow  interactions  occurring  as  the  outer  stream 
and  nozzle  exhaust  expand  into  the  lower  pressure  of  the 
base,  mix  with  the  base  gas,  and  recompress  to  the  down- 
stream pressure  are  not  \>rell  understood.  However,  it  is 
of  importance  to  the  designer  to  know  the  properties  of  the 
mixture  of  the  low-velocity  gases  recirculated  into  the 

base  region  as  they  are  at  a pressure  lower  than  the 

i 

ambient  stream,  and  they  contribute  significantly  to  the 
drag  of  the  vehicle.  Thus,  a detailed  knowledge  of  how  the 
base  pressure  and  temperature  are  affected  by  the  boundary 
layer  buildup  on  the  outer  surface  of  the  afterbody,  by  the 
flare  angle  of  the  nozzle,  by  the  constituents  of  the 
nozzle  gas,  and  by  the  injection  of  various  gases  into  the 
base  is  required  in  order  to  optimize  a design.  Those 
interested  in  the  structure  of  rocket  plumes  are  concerned 
with  the  effect  base  properties  have  on  the  Mach  disk 
location  and  the  destruction  of  the  coherent  inviscid 
plume.  And,  of  course,  the  engineer  in  planning  a test 
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Figure  1.  General  plume-stream  interaction 
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configuration  must  have  some  indication  of  what  to  expect 
from  a given  vehicle  design  at  various  simulated  flight 
conditions,  a method  of  confirming  anomalies  in  his  data, 
and  the  knowledge  to  interpret  and  extend  the  usefulness  of 
his  data. 

It  is  the  purpose,  then,  of  this  investigation  to 
present  a new  theoretical  model  capable  of  predicting  the 
base  properties  behind  blunt-based  bodies  with  supersonic 
nozzles  at  supersonic  speed  taking  into  account  the 
boundary  layer  existing  at  separation,  base  bleed,  and 
chemistry.  This  study,  however,  must  be  considered  a 
progress  report  since,  as  is  always  the  case,  limitations 
of  time  and  resources  have  left  several  avenues  of  investi- 
gation underdeveloped. 

Although  this  work  is  new  and  comprehensive,  it 
makes  use  of  well-established  principles  first  presented  by 
H.  H.  Korst  and  his  co-workers  [ 1] 1 and  extended  by  Bauer 
and  Fox  [2).  Thus,  it  is  perhaps  appropriate  to  focus  on 
the  primary  departures  from  earlier  concepts. 

One  of  the  primary  obstacles  to  a comprehensive 
base-flow  theory  is  the  recompression  process.  If  the 
classical  planar,  backward-facing  step  is  considered,  for 
example,  the  recompression  process  serves  to  compress  the 
high-speed  gas  in  the  mixing  layer,  from  what  is 


''“Numbers  in  brackets  refer  to  similarly  numbered 
references  in  the  Bibliography. 
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essentially  the  pressure  existing  within  the  base  region, 
to  the  free-stream  pressure.  If  the  flow  were  inviscid, 
the  pressure  would  be  adjusted  abruptly  by  an  inviscid 
shock.  However,  experiments  show  that  viscous  effects 
stretch  out  the  process  significantly.  The  streamline 
dividing  that  part  of  the  mixing  layer  which  experiences 
recompression  from  that  part  which  turns  back  toward  the 
base  must  be  identified  in  order  for  the  theory  to  have 
closure.  It  is  identified  by  the  pressure  to  which  it 
stagnates.  As  this  pressure  is  unknown,  earlier  workers 
made  up  empirically  determined  recompression  factors  to 
provide  the  required  information.  In  the  present  work,  a 
hypothesis  is  put  forward  which  serves  to  couple  the 
momentum  of  the  oncoming,  viscous  layer  to  the  observed 
recompression  profile  of  the  pressure  on  the  downstream 
wall.  This  hypothesis  is  a fundamental  refinement  to  the 
geometrical  recompression  theory  of  Bauer  and  Fox  [ 2] . 

When  considering  base  bleed  or  the  general  problem 
of  the  mixing  of  disparate  streams  into  the  base  region, 
the  distribution  of  the  species  through  the  mixing  zone 
presents  difficulties  for  existing  theories  when  the 
initial  boundary  layer  is  taken  into  account.  The  usual 
procedure  is  to  consider  the  species  profile  to  be 
identical  to  the  velocity  profile.  To  understand  that  this 
is  incorrect,  one  needs  only  to  consider  that  at  the  point 
of  separation  the  species  profile  is  that  of  a step 
function  while  the  velocity  profile  is  that  of  the  initial 
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boundary  layer.  Consequently,  in  the  present  theory,  a new 
profile  for  the  species  distribution  is  proposed  which 
serves  to  account  appropriately  for  this  distinction.  This 
profile  nicely  approaches  the  velocity  profile  in  the  limit 
of  fully  developed  flow. 

When  extending  the  theory  to  axially  symmetric 
flows,  the  mixing  theory  must  be  modified.  An  analysis  is 
presented  in  this  work  which  demonstrates  that  the  Mangier 
transformation  serves  to  confirm  the  mixing  theory  of  Bauer 
and  Fox  [2] . 

And  finally,  since  the  mechanisms  were  developed 
with  the  overall  theory,  it  was  a relatively  simple  matter 
to  incorporate  equilibrium  chemistry  into  the  analysis. 

Thus  the  use  of  the  so-called  flame-sheet  model  is  avoided. 

Although  .the  present  work  considers  only  turbulent 
flow,  the  extension  to  laminar  flow  is  straightforward. 

The  principal  change  would  be  in  the  analysis  of  the 
spreading  parameter. 

II.  Primary  Analytical  Models 

Even  with  the  advent  of  the  new  generation  of 
advanced  computers,  attacking  such  formidable  problems  as 
base  flow  with  bleed  and  chemical  effects  by  numerically 
solving  the  Navier-Stokes  equations  coupled  with  equi- 
librium chemistry  is  beyond  present  capabilities.  Conse- 
quently, the  approach  has  been  and  will  continue  to  be  that 
analagous  to  the  boundary- layer  problem;  that  is,  to  solve 


5 


A EDC-TR  -79*46 


the  viscid  and  inviscid  areas  by  different  methods  and  to 
couple  the  solutions  using  the  physical  constraints  at  the 
interfaces . 

There  are  two  primary  mixing  theories  extant  that 
may  be  used  in  the  base-flow  problem:  the  Crocco-Lees  [3] 
theory  and  the  component  model  of  Korst,  et  al . [1].  The 
Crocco-Lees  theory  and  the  Korst  theory  are  similar  only  in 
one  aspect:  they  both  account  for  the  strong  viscous  inter- 
actions that  prevail  in  flows  with  mixing. 

Crocco-Lees  Model 

The  basic  boundary- layer  equations  are  integrated 

across  the  mixing  layer.  There  results  a set  of  integro- 

differential  equations  which  are  ordinary  in  the  axial 

direction.  The  equations  are  linear  with  respect  to  the 

derivatives  at  any  axial  station.  Thus,  with  the  edge  con- 
\ 

ditions  known  and  the  transverse  profiles  given,  the  mixing 
layer  is  calculated  in  a marching  fashion  by  solving  the 
inviscid  field  simultaneously  to  provide  the  edge 
conditions.  Thus,  the  axial  pressure  gradients  are  part  of 
the  solution.  In  the  flow  of  a single  stream  over  a bluff- 
based  body,  the  elliptic  nature  of  the  base  flow  is 
exhibited  in  an  interesting  fashion:  if  the  guess  of  the 
initial  base  pressure  is  incorrect,  the  determinant  of  the 
derivatives  goes  to  zero  at  some  point  in  the  recompression 
region.  By  properly  adjusting  the  base  pressure,  the 
marching  of  the  solution  is  made  to  ptoceed  smoothly 
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through  this  point  with  the  determinant  of  the  augmented 
matrix  going  to  zero  simultaneously  with  that  of  the  coef- 
ficient matrix.  The  method  is  physically  realistic  and 
includes  the  details  of  the  actual  flow. 

The  main  drawbacks  to  this  method  which  foreclose 
its  use  in  the  present  analysis  are  twofold.  Of  first  con- 
sideration is  the  complexity  and  difficulty  of  imple- 
mentation of  the  method  in  any  but  the  simplest  of  problems; 
it  is  not  clear,  for  example,  when  there  are /two  supersonic 
streams  involving  two  converging  mixing  layers  as  well  as  a 
bleed  flow,  just  how  the  base  pressure  is  determined.  Two 
separate  singularities  would  have  to  be  coupled  and 
resolved  requiring,  at  least,  a detailed  analysis  of  the 
reversed  flow,  such  analysis  making  the  Navier-Stokes 
approach  a delight  by  comparison.  Secondly,  and  equally 
important,  is  the  finding  by  Peters  and  Phares  [4]  that  the 
method  breaks  down  when  the  species  concentration  equation 
is  coupled  into  the  system:  anomalous  singularities  occur. 
(The  problem  apparently  can  be  overcome  by  uncoupling  the 
species  equation  from  the  system  by  appropriate  modeling.) 

'i 

These  concerns  make  the  Crocco-Lees  approach  unattractive 
at  this  stage  of  its  development. 

Korst  Mixing  and  Base-Flow  Models 

The  basic  Korst  theory  is  so  uncomplicated  that  for 
simple,  fully  developed,  two-dimensional  base  flows  the 
computations  required  can  be  performed  on  a desk  calculator. 
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The  basic  simplicity  of  the  mixing  theory  along  with  the 
good  correlations  that  result  from  its  use  have,  made  .it  one 
of  the  most  widely  adopted  approaches  to  turbulent  mixing 
problems  [5].  A recent  example  of  its  use,  essentially 
unmodified,  is  afforded  by  the  work  of  Yaros  [6] , who  used 
it  to  account  for  mixing  between  two  coaxial  streams.  His 
results  agreed  well  with  observation. 

The  elemental  Korst  base-flow  theory  is  made  up  of 
three  distinct  analyses : 

1.  Inviscid  analysis  to  determine  the  overall  flow 
structure  and  pressure  field,  sometimes 
referred  to  as  the  "corresponding'’  inviscid 
flow. 

2.  Global  analysis  of  the  assumed  flow  behavior  in 
the  near  wake  leading  to  the  primary  system  of 
conservation  equations  for  the  base  properties. 

3.  Viscous  mixing  and  recompression  analysis  to 
determine  key  streamlines  in  the  mixing  region. 

The  simplicity  is  achieved  through  the  assumption 
of  isobaric  mixing.  This  allows  the  mixing  region  to  be 
completely  determined  at  any  point  along  the  inviscid 
boundary  as  a function  of  the  mixing  pressure  and  the 
distance  from  inception.  Thus  the  mixing  analysis  is 
uncoupled  from  the  inviscid  computations.  The  velocity 
profile  is  determined  by . s implif ication  of  the  boundary- 
layer  momentum  equation  to  a heat  conduction  form  that  has 
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an  error  function  solution.  The  only  coupling  to  the 
inviscid  flow  is  achieved  by  positioning  the  layer  with 
respect  to  the  inviscid  boundary  by  a momentum  balance. 

For  the  simple,  two-dimensional,  base-pressure 
problem  (flow  over  a backward- facing  step),  the  Korst 
theory  assumes  that  the  dividing  streamline,  which 
discriminates  between  that  part  of  the  flow  which  goes 
through  recompression  and  that  which  recirculates  in  the 
base  region,  isentropically  stagnates  to  a pressure  equal 
to  the  recompressed  pressure.  This  knowledge  allows  a 
global  mass  conservation  equation  for  the  recirculating 
flow  to  be  evaluated.  The  pressure  is  adjusted  until  a 
dividing  streamline  is  determined  that  allows  the  global 
mass  balance  equation  to  be  satisfied.  At  first  glance, 
this  appears  to  be  an  oversimplification  of  the  problem; 
however,  it  models  the  major  elements  of  the  flow  and  corre- 
lates well  with  experiment. 

Further  investigation  led  to  the  conclusion  that 
the  remodeling  of  elements  of  the  method  would  make  it  more 
general  and  thus  applicable  to  the  complex  flows  under  con- 
sideration. Therefore,  the  Korst  mixing  and  base-pressure 
theory,  as  developed  and  extended  in  the  following  chapters, 
is  the  method  of  choice. 
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III.  Other  Approaches 

Over  the  years,  the  Korst  and  Crocco-Lees  theories 
have  stimulated  similar  approaches  in  an  effort  to  find 
better  agreement  with  observation.  Several  are  reviewed. 

Chow's  Method 

Chow  [7]  developed  a base-pressure  theory  which 
uses  the  Korst  constant-pressure  mixing  up  to  the  point 
where  recompression  begins;  whereupon,  the  mixing  layer  is 
split  at  the  dividing  streamline,  the  outer  part  exposed 
to  a pressure  gradient.  The  high-speed  layer  goes  through 
recompression,  while  the  low-speed  layer  is  obverted  and 
assumed  to  recirculate.  The  shear  force  along  the  dividing 
streamline  must  be  equal  at  once  to  the  gradient  of  the 
momentums  in  each  layer.  By  modeling  the  shear  stress,  the 
gradients  are  known,  and  an  integration  in  the  direction  of 
principal  flowr  can  be  performed  stepwise  to  the  stagnation 
point,  where  the  velocity  of  the  dividing  streamline  is  set 
to  zero  and  the  pressures  assume  particular  values  derived 
from  momentum  considerations  on  the  stagnation  surface.  By 
assuming  a base  pressure  and  the  distance  from  separation 
at  which  the  recompression  begins,  the  integration  is  able 
to  proceed  to  stagnation  where  the  dependent  variables  must 
match  the  boundary  conditions  on  the  stagnation  surface. 

If  they  do  not,  a new  guess  of  pressure  and  length  is  made 
and  an  iteration  ensues  until  the  integration  predicts  con- 
ditions at  stagnation. 
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This  analysis,  while  containing  \rery  realistic 
assumptions,  is  fairly  complex  and  gives  only  fair  pre- 
dictions for  thin  boundary  layers  and  poor  results  for  very 
thick  boundary  layers. 

Chow  and  Spring's  Modification 

This  theory  [8]  is  essentially  a modification  of 
the  preceding  theory.  The  constant -pressure  mixing  up  to 
the  onset  of  recompress ion  is  replaced  with  a refinement 
which  includes  a streamtube  expansion  of  the  initial 
boundary  layer  at  the  corner  separation  point.  A linear 
velocity  profile  is  then  assumed  that  has  the  slope  of  an 
error -function  profile  at  its  half -velocity  point.  The 
method  of  solution  is  quite  similar  to  the  previous  theory. 
The  predictions  of  the  theory  are  quite  good  for  thicker 
boundary  layers;  however,  here  again  the  method  is  quite 
cumbersome . 

Tanner's  Method 

Tanner's  [9]  approach  is  more  analogous  to  the 
Korst  theory  than  any  other  and  is  therefore  attractive  in 
its  simplicity.  His  principal  assumption  is  that  the  pro- 
duction of  entropy  of  the  recompression  shock  in  an 
inviscid  flow  is  identical  to  the  entropy  production  when 
the  reattached  shear  layer  displaces  the  shock.  That  is, 
the  lost  entropy  production  of  the  displaced  shock  is  made 
up  by  the  entropy  production  of  the  shear  layer . Tanner 
determines  empirically  the  ratio  of  the  distance  from  the 
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displaced  inviscid  shock  and  the  thickness  of  the  recom- 
pressing  shear  layer.  He  also  assumes  the  velocity  at  the 
stagnation  surface  after  the  flow  has  reaccelerated.  Base  . 
pressures  are  guessed  in  an  iterative  fashion  until  the 
condition  of  entropy  equality  is  satisfied.  Although 
Tanner  does  not  take  into  account  the  initial  boundary 
layer,  his  agreement  with  his  own  data  is  quite  good;  in 
fact,  the  theory  of  this  investigation  agrees  with  that  of 
Tanner  for  a particular  boundary- layer  momentum  thickness 
(see  Chapter  IV,  Section  II).  The  principal  objection, 
however,  is  that  it  is  not  clear  just  how  bleed  flow  could 
be  incorporated  into  the  theory,  as  Tanner  does  not  employ 
any  mass -rate  accounting. 

Extensions  of  the  Basic  Korst  Model 

Nash  [10],  Addy  [11],  and  Page,  et  al . [12]  are 
examples  of  the  second  generation  of  Korst  theorists.  They 
all  resort  to  an  empirically  determined  recompression 
factor  to  make  their  solutions  agree  with  various  data. 
Thus,  as  Peters  and  Phares  [4]  point  out,  ".  . . these 
models  really  are  used  to  correlate  experimental  base- 
pressure  measurements  rather  than  to  predict  them,"  as  no 
universal  factor  has  been  developed. 

Refinements  of  the  Crocco-Lees  Approach 

The  works  of  Peters  and  Phares  [4] , Alber  and 
Lees  [13] , and  Strahle  and  Mehta  [14]  present  extensive 
refinements  to  the  basic  Crocco-Lees  approach.  Although 
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each  differs  in  significant,  detail,  such  as  the  determi- 
nation of  the  turbulence  - transport  model  and  the  treatment 
of  the  inviscid  stream,  they  all  exhibit  to  varying  degrees 
the  difficulties  of  the  basic  theory  as  pointed  out  in 
Section  II.  They  all,  for  the  most  part,  give  good  results 
for  simple  flows.  The  extension  to  two  converging  mixing 
layers  is  not  straightforward. 


\ 
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CHAPTER  II 

THEORY 

I . Basic  Considerations 

When  a turbulent  stream  separates  at  the  bluff 
trailing  edge  of  a body,  the  free  tangential  shearing  sur- 
face, being  unstable,  breaks  into  eddys  which  transport 
matter  across  the  surface  creating  a mixing  zone.  . The 
mixing  zone  is  an  identifiable  region  across  which  the 
properties  of  the  stream  and  base  gases  are  distributed. 
Time-averaged  measurements  show  these  distributions  as 
coherent  profiles  of  velocity,  temperature,  and  species 
concentration.  (See  Hussain  and  Zedan  [IS]  and  Rhudy  and 
Magnan  [16].)  After -turning  toward  the  lower  pressure  of 
the  base,  the  stream  must  recompress  and  return  to  its 
original  direction.  Part  of  the  mixing  layer  will  have 
insufficient  energy  to  overcome  the  adverse  pressure  change, 
and  it  will  remain  in  the  recirculating  fluid  of  the  base. 

If  a pure  gas,  say  hydrogen,  were  bled  into  the  base,  the 
base  gas  would  be  a mixture  of  hydrogen  and  the  free-streain  . 
air.  In  steady  flow,  the  base  gas  can  be  considered  as  a 
separate,  quiescent,  homogeneous  mixture.  The  mixing  layer 
would  be  a distribution,  then,  from  pure  base  mixture  to 
pure  free-stream  air  at  the  high-speed  edge.  At  an  inter- 
mediate point  within  the  layer,  the  local  mixture  would  be, 
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say,  40  percent  base  mixture  and  60  percent  free-stream  air. 
The  base  mixture,  of  course,  would  contain  some  free-stream 
air . 

Reynol ds - averaged  flow  is  assumed  throughout  the 
following  analytical  development. 

II.  Key  Streamlines 

Two  streamlines  which  are  important  to  the  evalu- 
ation of  the  global  conservation  equations  are  the 
stagnating  streamline,  yg , and  the  dividing  streamline,  y^. 

The  stagnating  streamline  is  that  streamline  which 
separates  that  part  of  the  flow  which  continues  through 
recompression  (the  high-speed  part)  from  that  which  turns 
back  into  the  base  region  (the  low-speed  part) . It  is 
determined  by  the  recompression  theory  of  this  chapter  in 
Section  VIII. 

The  dividing  streamline  is  defined  simply  for  mass 
accounting  purposes.  It  is  defined  by  the  following 
relation: 

yU  yu 

f p_u  dy  ] = f pudy  ] . (1) 

o x=o  y D x=£ 

That  is,  the  dividing  streamline  is  that  streamline  outside 
of  which  the  mass  of  the  separating  flow  at  the  corner  is 
completely  accounted  for. 
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III.  Open  Wake 

The  open  wake  was  hypothesized  by  Korst,  et  al . [1], 

It  serves  to  account  for  the  fact  that  a net  mass  of  fluid 
is  being  pumped  into  and  out  of  the  base  region  by  the 
high-speed  streams.  With  reference  to  Fig.  2,  if  the 
dividing  streamline  were  found  on  the  low-speed  side  of  the 
stagnating  streamline,  then  part  of  the  mass  of  the  stream 
before  separation  is  counted  as  going  into  the  base.  The 
converse  situation  would  have  part  of  the  base  flow 
eventually  being  ushered  downstream  through  recompression. 

In  the  limiting  case  of  symmetrical,  two-dimensional  flow 
without  base  bleed,  the  dividing  streamline  and  the 
stagnating  streamline  coincide. 

IV.  Global  Conservation 

The  global  mass  balance  may  now  be  formulated.  In 
other  words,  the  mass  rate  of  the  base  flow  entrained  into 
the  mixing  layer  must  equal  the  mass  rate  of  that  which  is 
turned  back  into  the  base  region,  plus  any  bleed. 

Consider  the  control  volume,  A-B-C-D,  of  Fig.  2 
which  encloses  the  base  region.  The  mass  flowing  inside 
the  dividing  streamlines  accounts  for  all  the  mass  of  the 
base  region.  The  mass  crossing  A-B.  and  C-D  is  the 
entrained  base  mass;  so  that 
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Figure  2.  Schematic  of  flow  behavior  in  near-wake  region, 
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yD  yD 

B + nip  D = y pudy  J + f pudy  J = ‘ (2) 

yL  S1  y L S2 

The  mass  entering  B-C  is  accounted  for  as  the  total  mass 
flow  turned  back  inside  the  stagnating  streamlines,  or 

yS  yS 

mr  n = / pudv  ] + / pudy  ] . (3) 

yL  31  yL  32 


With  bleed  flow,  the  mass  conservation  equation  becomes 


yD  yD  yS  yS 

/ pudy  ] + / pudy  ] = mF  + / pudy  ] + / pudy  ] (4) 

yL  SI  yL  S2  yL  SI  yL  S2 


Similarly,  the  energy  equation  is 


yD  yD  yS  yS 

Hi>{/  pudy]  + / pudy]  }=  IT1PHP+  / puHdy  ] + / puHdy  ] , (5) 

yL  SI  yL  S2  yL  SI  yL  S2 

where  the  mass  through  the  base  region  is  assumed,  as 
mentioned,  to  be  a homogeneous  mixture  of  low-velocity  gas. 
As  the  base-region  gas  and  the  bleed  are  assumed  to  have 
negligible  velocities,  there  are  sufficient  conservation 
relations  for  finding  two  thermodynamic  properties  of  the 
base  such  as  pressure  and  enthalpy.  However,  at  this 
point,  the  limits  of  integration  are  unknown,  as  are  the 
relationships  between  the  integrands  and  the  variables  of 
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integration.  The  appropriate  relations  will  be  developed 
in  the  following. 

V.  Mixing  Analysis 

Here  the  mixing  analysis  of  Korst  is  followed  as 
detailed  by  Chapman  and  Korst  [17],  In  particular,  the 
general  case  with  initial  boundary  layer  is  considered. 

Consider  the  equation  of  motion  of  a two- 
dimensional,  constant  pressure  boundary- layer  flow. 


Pu 


9u 

3x 


+ 


pv 


= A f £0  — 1 

By  3y  3y 


(6) 


In  turbulent  flow,  u and  v are  the  time -averaged  velocity 
components,  and  t is  interpreted  as  an  apparent  kinematic 
viscosity.  The  equation  of  motion  is  reduced  using  order 
of  magnitude  arguments  and  the  assumptions  that  the 
apparent  viscosity  is  independent  of  y and  the  flow 
incompressible;  that  is, 

2 

3u  _ e 5 u j-y-j 

9X  U gy2 

By  assuming 

e = e0FU)  , 

Eq.  (7)  may  be  cast  in  the  form 

? 

3u  = 3_u  . (9) 

3?  9y2 
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The  solution  to  the  heat  conduction  equation  is  well  known 
(see  [18] ) , It  is  : 

00  . (y-g)2 

u(y,C)  = / e f0(3)dS  , (10) 

where  £q(3)  is  the  initial  condition.  From  Fig.  2,  it  is 
apparent  that  the  initial  profile  is  zero  from  minus 
infinity  to  zero.  From  zero  to  the  boundary  - layer 
thickness,  5,  it  is  the  boundary- layer  profile.  From 
6 to  +°°,  the  profile  is  that  of  the  expanded  free  stream. 
Applying  this  initial  profile, 

fi  _ (y-e)2 

u(y, O = r-  ■[  1 + erf  [ ] } + — - — / e u (3)d0  (11) 

2 2ft.  2ffZ  Q S 

is  obtained. 

It  is  well  known  that  when  the  mixing  region  is 
fully  developed,  that  is,  when  S/x  is  very  small,  that  the 
velocity  profile  is  well  approximated  by 

u = ?[l  + erf  (H£)  ] , (12) 

where  a is  a constant  usually  referred  to  as  the  spreading 
parameter,  which  is  experimentally  determined.  Using  this 
fact,  the  velocity  profile  may  be  written  as 


i 
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u(n)  = 7 [ 1 + erf  ( n-n  ) ] +■  S e 

2 p sa 

0 


1.  '-0-  m-n „C?)I 


p 4 \cjf34(?3.  a3, 


where 


and 


n 


= £1 


X 


(14) 


This  implies  that  the  kinematic  viscosity,  e,  is  a linear 
function  of  x. 


VI.  Determination  of  y 

' m 

The  limits  of  integration  y^  and  yL  as  used  in 
Eq.  (1)  are  functions  of  distance  from  the  inception  of 
mixing.  The  upper  limit  is  taken  as  that  point  where  the 
velocity  ratio,  u/U,  is  0.99998895.  The  lower  limit  is 
taken  where  u/U  is  0.00001105.  Thus,  the  actual  y value 
changes  as  the  mixing  zone  spreads  into  the  inviscid  flow. 
At  a particular  distance  from  the  inception  of  mixing,  if 
the  particular  inviscid  streamline  at  y^  were  traced  back 
to  the  beginning,  the  momentum  of  the  mixing  zone  should 
account  for  the  momentum  of  the  initial  profile  out  to  that 
streamline.  Because  of  simplifications  of  the  theory,  how- 
ever, the  momentums  will  not  check.  So  as  to  balance  the 
momentum,  the  upper  limits  on  the  integrals  are  extended  by 
V . which  is  not  a coordinate,  but  an  increment;  so  that 
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s 


2dy  ] 
x=o 


Vy„ 

/ purely  ] 

yL  ■ x-i 


Korst  and  his  co-workers  also  consider  the  profile  to  have 
shifted  by  a like  amount  so  that  the  inviscid  boundary  no 
longer  is  located  at  the  half-velocity  point  of  the  fully 
developed  profile,  but  at  y = y . As  this  shift  causes  the 
profiles  to  align  better  with  experimental  profiles,  this 
localization  with  respect  to  the  inviscid  boundary  is 
adopted  here.  The  transverse  coordinate  is  always  fixed  at 
the  half-velocity  point  of  the  fully  developed  profile 
wherever  the  profile  may  fall.  Thus,  y is  completely 
determined  by  Eq . (15). 


VII,  Constitutive  Equations 


Mass  Fractions  of  the  Streams  in  the  Base 

To  be  able  to  determine  the  base  mixture,  the  mass 
fractions  of  the  primary  streams  and  the  bleed  gas  in  the 
base  must  be  determined.  As  the  base  mixture  is  determined 
by  the  fact  that  distributed  mixtures  (the  mixing  layers) 
are  being  turned  back  into  the  base,  the  functional  form  of 
the  species  distribution  must  be  determined  before  pro- 
ceeding. 

If  the  assumption  of  no  initial  boundary  layer  is 
made,  as  in  [19],  then  by  assuming  unity  turbulent  Prandtl 
and  Lewis  numbers,  the  species  and  total-enthalpy- distri- 
bution profiles  can  be  taken  as  the  same  as  the  velocity 
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profile.  However,  in  the  situation  being  investigated  here, 
the  initial  species  profile  is  a step  function;  whereas  the 
velocity  profile  is  that  of  the  boundary  layer.  By 
assuming  a species  distribution  function  as 

k = | [ 1 + erf(Cmfn)  ] , (16) 


with  the  shape  parameter,  Cmf>  to  be  determined,  the  dis- 
parity between  the  profiles  is  accommodated.  The  species 
profile,  which  has  larger  gradients,  is  always  approaching 
the  velocity  profile  as  the  mixing  develops.  (This  is 
borne  out  in  the  computations  which  calculate  C ^ to  be 
always  greater  than  unity  except  in  the  fully  developed 
limit  where  it  is,  of  course,  unity.) 

The  C ^ parameter  is  determined  by  a mass  balance, 
as 

yU  Vym 

/ p u dy  ] = / kpudy  ] . (17) 

o s x^o  yL  x=£ 


Thus,  with  Eq . (15)  as 

yu  , 

f PcUc“dX  1 
o x=o 


>U+>m  ? 

/ puzdy  ] 
yL  x-t 


CIS) 


there  are  two  equations  for  the  two  unknowns,  Cm^  and  ym . 
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Then  the  mass  fraction  of  Stream  1 gas  in  the  base 
mixture  is 


rS  Yg  Yg 

f kpudy  ] + mf,  { / (1-k)  pudy  ] + / (1-k)  pudy  ] } 

yL  S1  >l  S1  yL 

>’s  ys 

nip.  + / pudy  ] + / pudy  ] 

yL  SI  yL  S2 


or 


(19) 


yS 

/ kpudy  ] 


ys  ys 

rip  + / kpudy  ] + f kpudy  ] 

yL  SI  yL  S2 


Also  for  Stream  2 , 


mf2 


ys 

/ kpudy  ] 

yL  52 
ys  ys 

m-n  + / kpudy  ] + I kpudy  ] 

yL  S1  yL  52 


(20) 


(21) 


and  for  the  bleed  gas, 

I 


mfE 


1 - mf^  - mf2  . 


Thus , having  the  mass  fractions  of  the  respective 
stream  in  the  base,  sufficient  information  is  available  to 
determine  the  base-gas  mixture  and  its  properties. 
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Chemical  Equilibrium  Considerations 

In  any  compressible  flow  analysis,  a thermal 
equation  of  state  is  required  to  close  the  system  of 
.equations,  such  as 

P = P(p,h)  . (22) 

A caloric  equation  of  state  is  also  required  since 
empirically  determined  properties  are  expressed  as 
functions  of  temperature,  so  that 

t = t(p,h)  , - (23) 

The  requirement  of  the  equilibrium  chemistry  analysis,  then, 
is  that  it  supply  the  required  functional  relations. 

The  mixture  is  assumed  to  be  composed  of  thermally 
perfect  gases;  so  that  the  thermal  equation  of  state  is 

l ci«i 

P = £ ] . (24) 

R 

The  caloric  equation  is  of  the  form 


i 


which  assumes  that  each  species  behaves  ideally.  The  zero 
superscript  denotes  a known  evaluation  of  the  variable  at 
reference  (zero  temperature)  conditions,  To  evaluate  these 
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relations 
tion  of  t 
present . 


requi res 
emperature 


and 


for  each  species  present  as  a func- 
the  mole  fraction  of  each  species 


Reference  [20]  provides  a functional  form  for  the 
various  properties  of  210  substances  to  6000°K.  For 
example,  the  specific  heat  at  constant  pressure  is  given  as 


c 

Pi  7 7 4 

— = a1  + a2t  + a_t“  + a^t  + ast  , (26) 


But  to  make  the  problem  simpler,  this  reference  also  pro- 
vides, directly,  h^  = lu(t)  as  well  as  entropy  and  the 
Gibbs  function  per  mole.  All  that  is  further  required, 
then,  are  the  mole  fractions,  C,  of  the  species  present. 

In  general,  in  a given  mixture  of  m atoms  and 
n species,  the  following  relations  hold.  For  example,  see 
[21].  For  a given  mixture  there  are  n-m  equilibria 

y a..X.  i 3. .X.  i = 1,  n-m  , (27) 

h ii  j i]  J 


with  the  corresponding  mass -action  relations  as 


(3- • -a. .) 


Ki  = 11  1J  1]  = exP  I-1 


-I  CPy-ayiPy  i=l,  n-m 

— — ] j=l > n , (28) 


Rt 
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and  m atom  conservation  equations 


j akjpkj 


Jf 


k=l , m 

j=l  , n 


(29] 


where  is  the  fictitious  partial  pressure  of  the  kth 
atom. 

The  Gibbs  free  energy  per  mole  of  individual 
species  is  given  by  [20]  in  functional  form  as  mentioned. 
These  n equations  are  sufficient  to  solve  for  the  n partial 
pressures  of  species  present  and,  along  with  the  known 
pressure  of  the  mixture,  to  determine  the  mole  fractions  of 
the  species.  It  is  possible,  then,  for  a given  mixture  of.  • 
atoms  to  determine  any  state  variable  required  as  a 
function  of  pressure  and  enthalpy.  Appendix  A show's  by 
example  the  method  of  determining  the  mole  fractions. 

Flame  Sheet 

At  this  point,  it  is  necessary  to  examine  the 
flame-sheet  concept  as  implemented  by  Davis  [19]  . This 
concept  is  attributed  to  Abramovich  [22],  and  Libby  [23]. 
The  primary  assumption  of  the  model  is  that  combustion 
occurs  along  an  infinitely  thin  layer  within  the  mixing 
region.  At  this  point,  the  concentrations  of  fuel  and 
oxidizer  are  taken  to  be  zero,  but  the  amounts  consumed  in 
the  reaction  are,  in  stoichiometric  proportion.  The  con- 
centration of  fuel  on  the  lean  side  of  the  flame  is  zero  as 
is  the  oxidizer  on  the  rich  side.  In  Davis'  case,  he 
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considers  only  the  hydrogen/oxygen  reaction  of  his  experi- 
ment. By  assuming  a Crocco  relation  between  the  velocity 
profile  and  the  atomic  species  distribution,  he  is  able  to 
determine  the  stoichiometric  point  in  the  mixing  layer, 
then  distribute  the  products  of  combustion  on  either  side 
of  the  flame  front  using  the  Crocco  relation  and  the  con- 
centration assumption  on  the  fuel  and  oxidizer  as  stated 
above.  Although  no  motivation  is  given  by  Davis  for  this 
approach,  it  is  quite  apparent  that  the  need  to  solve  the 
mass-action  equations  and  atom- conservation  equations  for 
the  equilibrium  concentrations  of  each  species  has  been 
obviated.  The  principal  objection,  though,  to  Davis'  ap- 
proach is  that  he  disallows  oxygen  to  be  present  in  the 
base  mixture.  This  would  bring  into  question  Davis'  analy- 
sis of  the  cases  with  little  hydrogen  bleed,  which  have  base 
mixtures  on  the  lean  side  of  the  stoichiometric  point. 

In  conclusion,  it  can  be  said  that  Davis'  chemistry 
model  is  a relatively  good  approximation  for  mixing  layers 
between  streams  of  pure  fuel  and  oxidizer.  In  the  base- 
flow  problem,  these  conditions  do  not  obtain  because  of 
the  recirculation  into  the  base  of  significant  portions  of 
the  oxidizer  stream. 

VIII.  Recompression  Analysis 
General  Considerations 

Korst's  original  analysis  did  not  contain  a recom- 
pression model  as  such;  rather,  he  simply  assumed  that  the 
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stagnating  or  discriminating  stream  isentropical ly 
stagnated  at  the  static  pressure  downstream  of  the  recom- 
pression shock.  Experiments  showed  this  to  be  incorrect. 
To  account  for  this  fact,  others  such  as  Nash  [10],  and 
Addy  [11]  devised  empirically  determined  recompression 
factors.  These  were  unsatisfactory  as  they  could  only 
correlate  existing  data  rather  than  predict. 

Bauer  and  Fox  [2],  however,  devised  an  analytical 
model  which  obviated  the  requirement  for  a recompression 
factor.  Their  approach,  however,  overpredicted  the  base 
pressure  for  thick  boundary  layers.  This  w'as  because  they 
used  essentially  an  inviscid  analysis  to  determine  the 
extent  of  the  recompression  region. 

In  the  present  work,  a novel  approach  is  used  for 
modeling  the  recompression  process. 


Approach 

Determination  of  the  location  of  the  stagnating 
streamline,  ygt  is  tied  to  the  recompression  process  where 
the  two  streams  merge,  turn,  and  recompress  to  the  ambient 
pressure.  Figure  3 isolates  one  of  the  impinging  streams. 

The  base  flow  is  assumed  to  turn  back  through  an 
area  equal  to  the  approaching  area  between  yg  and  yL;  thus 
the  beginning  of  recompression,  X^  is  related  to  yg 
geometrically.  Reference  [2]  showed  the  function,^ 


P ~ PB 
Pr  -pB 


sin 


X - X, 

T - 1 


It  ( 


X2~X1 


)] 


(30) 
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to  be  a good  fit  of  recompression  data.  The  downstream 
pressure,  p , is  that  resulting  from  the  inviscid  calcu- 
lation of  shock  and  slipline.  The  actual  recompression 
process  on  the  slipline  is  stretched  out  as  indicated  by 
Eq.  (30).  The  end  of  the  recompression  process,  X_  , is 
approximately  determined  from  simple  momentum  consider- 
ations . 

It  is  assumed  that  the  force  due  to  the  pressure 
rise  on  the  slipline  is  equal  to  the  transverse  momentum 
given  up  by  the  oncoming  viscous  layer,  with  all  other 
effects  being  of  negligible  importance.  This  leads  to  the 
relation 


/ pdX  = sin  y / pu^dy  ] + ( xu  ' X1  ’ 

X1  XL  X=* 


where 

xU*Xm+  (yU  'xm)  1 sinY  ‘ 


(31) 


(32) 


For  the  backward- facing  step  problem,  Peters  and 
Phares  [4]  show  that  the  total  pressure  on  the  dividing 
streamline  changes  little  up  to  stagnation,  thus  the  stag 
nation  process  can  be  approximated  as  isentropic.  The 
stagnating  streamline  is  determined  by  finding  that 
stagnation  pressure  in  the  oncoming  mixing  layer  which  is 
equal  to  the  static  pressure  of  the  recompression  region  at 
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the  stagnation  point,  Xg.  The  recompression  pressure  at 
stagnation,  pg,  is 


also 


and 


Y - Y 

■ 2 r 7T  , S 1 % , 

Sln  [TcxrofT)J 


2 

hs  = Hs  - ug/2  , 
s = s ( PB , hg  3 , 


S — S ( Pg,  Hg  ) , 


(33) 


(34) 


where  Ug  and  Hg 
profiles,  and  y 
(Fig.  3).  The 
determining  the 
pressure . 


are  functions  of  yg  through  the  assumed 
s and  Xg  are  related  by  the  geometry 
above  relations,  then,  are  sufficient  for 
stagnating  streamline  as  a function  of  base 


IX.  Spreading  Parameter  for  Turbulent  Mixing 


General 

Perhaps  one  of  the  most  important  aspects  of  any 
mixing  analysis  is  the  determination  of  the  spreading  or 
mixing  parameter,  o.  As  Davis  [19]  points  out,  the  experi- 
mental data  are  exiguous  for  jet  mixing  with  combustion. 

For  isoenergetic  mixing,  the  most  widely  used  relation  is 
that  due  to  Korst  and  Tripp  [24], 
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a = 0Q  + 2.758  Mu  , (35) 

where  oq  is  the  incompressible  spreading  parameter. 
Channapragada  [25]  presented  a relation  using  a mixing 
length  argument  which  attempts  to  account  for  density  vari- 
ation as 


P P - 1 

a = a [ B ( U0+  B ) ] , (36) 

O Py 

where  B is  a correlation  factor  that  has  values  of  from  0.5 
to  0.25  as  Mach  number  increases  from  zero  to  above  3.0. 
Abramovich  presented  a relation  that  reduces  to 

a = 2aQ  ( 1 + P B /P jj  )~ 1 (37) 

for  mixing  with  a quiescent  base  gas. 

Prandtl’s  Mixing  Length  Model 

The  most  satisfactory  model  for  this  analysis,  how- 
ever, is  that  of  Bauer  [26]  who  ties  the  Prandtl  mixing- 
length  hypothesis  directly  to  the  integral  approach  of  the 
Korst  theory. 

Considering  the  fully  developed  profile  of  Fig.  4 
and  using  a control  volume  approach,  it  is  apparent  that 

v 

£.  ' D 

j'TydJl  = f pu2dy  ' (38) 

0 yL 

from  a momentum  balance.  Any  viable  shear  stress  model 
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could  be  used  for  r^.  However,  in  this  approach,  the 
Prandtl  mixing-length  model  is  used  as 

Tp  = pD^  2 ( du/ dy  ) ^ . (39) 

It  is  assumed  that  the  mixing  length  is  proportional  to  the 
local  mixing  zone  width: 

i = Zr\jjkl/c  (40) 

or 

td  = 4njk2  pD  ( du/dp  )2  • (41) 

Substituting  into  Eq.  (38)  and  performing  the  integration, 

"U  2 

O = / — 4'2dn  / I c,  ( ^-  ) ] (42) 

p k dn  D 

nL  u 

is  obtained.  The  constant,  Cj,  is  determined  by  [26]  as 
C .5085/ao. 

A comparison  of  the  various  theories  is  presented 
in  Fig.  5 with  cq  as  12.0.  As  is  pointed  out  by  Birch  and 
Eggers  [27],  the  data  follow  two  distinct  curves.  The  main 
difference  in  the  data  is  Reynolds  number:  the  lower  data 
curve  having,  in  general,  the  lower  Reynolds  number.  The 
best  model  for  the  present  analysis  is,  as  mentioned,  the 
lowest  curve  [26].  It  changes  the  least  with  Mach  number. 
It  seems  that  base  flows  in  particular  require  the  flat 
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Figure  5. 


Comparison  of  the  various  spreading  parameter 
models  , 
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trajectory.  It  is  as  though  the  entrained  flow  from  the 
base  retained  a memory  of  its  turbulent  past,  which  serves 
to  increase  mixing  against  the  opposite  trend  of  the  Mach 
number  effect  for  the  mixing  of  flows  which  have  not  been 
recirculated . 

It  must  also  be  pointed  out  that  the  reference 
value  of  12.0  is  for  developed  flows.  Most  base  flows  are 
not  fully  developed.  Hussain  and  Zedan  [15]  indicate  that 
the  history  of  the  initial  boundary  layer  has  a significant 
effect  on  the  developing  a . Their  experiments  show  that, 
for  a boundary  layer  that  has  been  tripped,  the  mixing 
parameter  takes  between  one  and  two  nozzle  radii  to  reach  a 
value  of  nine.  Consequently,  was  taken  as  nine  in  all 
flows  except  those  involving  bleed  with  burning  where  a 
value  of  twelve  was  used  because  of  convergence  problems 
with  the  value  of  nine. 

Using  a of  nine  with  the  mixing- length  theory 
gives  values  of  the  computed  mixing  parameter  that  range  in 
value  from  9.0  to  about  13.5  for  the  isoenergetic  problems. 
This  is  consistent  with  the  constant  value  of  12.0  which  is 
used  by  others  such  as  Strahle  and  Mehta  [14],  and  Chow  and 
Spring  [8].  In  fact,  a constant  value  of  12.0  would  work 
quite  well  for  the  present  work. 

X.  Correction  for  Axially  Symmetric  Effects 

The  theory  as  developed  is  for  two-dimensional 
flows.  Extending  the  Korst  analysis  with  initial  boundary 
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layer  to  axially  symmetric  flow  is  not  obvious.  In  par- 
ticular, the  procedure  for  determining  the  profile  shift  is 
not  developed  in  the  literature.  A suggested  approach  is 
developed  in  Appendix  B;  but,  at  this  point,  it  appears 
costly  in  computational  time. 

However,  as1  the  two-dimensional  theory  is  an 
approximation  and  gives  good  results,  it  is  reasonable  to 
attempt  to  approximate  the  axially  symmetric  effects.  The 
first  approach  attempted  was  that  used  in  [2], 

Bauer  [28]  showed  that  for  a zero  initial  boundary 
layer,  under  circumstances  that  prevail  in  many  of  the 
flows  of  interest,  the  two-dimensional  mixing  analysis  is  a 
good  approximation  for  axially  symmetric  flows  if 


oRj 

— „ > 10 

% cos  3 


C43) 


This  condition  holds  for  the  flows  under  consideration  in 
this  study.  I t is  reasonable  to  assume  that  the  inclusion 
of  the  initial  boundary  layer  would  not  alter  this  con- 
dition. 

Reference  [2]  also  used  a new  spreading  rate 
parameter  which  was  adjusted  to  account  for  the  difference 
in  spreading  between  axially  symmetric  mixing  and  two- 
dimensional  mixing.  This  is  derived  in  [2]  as 


2Ri 

°k~a  ^ R + R-  J 


(44) 
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where  the  bar  serves  to  distinguish  the  two-dimensional 
variables  for  this  discussion.  Reference  [2]  goes  further 
and  assumes  that  the  variation  in  radius  across  the  mixing 
layer  is  small  and  can  be  removed  from  beneath  the  integral 
sign.  This  has  the  effect  of  making  the  global  conser- 
vation equations  identical  to  those  in  two-dimensional  flow 
when  no  bleed  is  present.  The  theory  is  summarized  in  the 
following . 

1.  The  inviscid  field  is  determined  from  the 
axially  symmetric  method  of  characteristics  . 

2.  The  mixing  theory  and  the  recompression 
process  are  treated  as  two-dimensional  except 
as  modified  by  a new  a. 

3.  The  global  conservation  equations  are  two- 
dimensional  except  as  modified  by  the  new  a. 

While  the  assumption  of  essentially  two-dimensional 
mixing  may  be  justified  by  the  nature  of  the  mixing  and  the 
fact  that  it  is  removed  from  the  centerline,  it  is  not  as 
clear  that  the  mixing  layer  may  be  considered  thin  as 
regards  the  recompression  process  and  the  global  conser- 
vation . 

Consequently,  a second  approximation  was  imple- 
mented whereupon  the  global  conservation  relations  and  the 
recompression  relations  were  replaced  by  their  axially 
symmetric  counterparts.  Typically,  for  mass  conservation, 
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'D  /D  _ yS  JS 

2tt  / puRdy  ] + Zn  f puRdy  ] = m + 2tt  / puRdy  ] + 2tt  / puRdy  ] , (45) 

V ^1  V Q?  ^ ar  Cl  x t CO 


SI 


S2 


yL 


SI 


yL 


S2 


and  for  the  recompression  momentum  balance, 


yU+ym 


2tt  / pRdX=  2tt  sin  y / pu2Rdy  + it  Pg  ( + R^)  (X^  - X^) 

X1  y\. 


(46) 


A third  approximation  is  considered  which  is  based 
on  the  familiar  Mangier  transformation.  Consider  first 
that  the  entire  two-dimensional  approximation  is  based  on  a 
solution  of  the  simplified  boundary -layer  equations.  For 
laminar  flow,  the  well-known  Mangier  transformation  (see 
Tetervin  [29])  transforms  the  axisymmetric  boundary-layer 
equations  into  the  two-dimensional  boundary-layer  equations 
exactly.  For  a free  boundary,  the  transformation  is 


and 


y 

P u(x) 
p uU(x) 

pv 


A R.  M 
■ f C ) dx  , 

R. 

■f'  - 

: pu(x)  , 

PyU  (X)  , 


R. 

1 


r - dR, 

(p  v)  - p u , 


(47) 
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where  R_^  is  the  radius  of  the  inviscid  boundary  and  L is  a 
reference  length.  As  the  two- dimensional  analysis  gives 
good  results,  the  question  is:  would  this  simple  trans- 
formation give  a good  axisymmetric  approximation? 

Reference  [29]  points  out  that  the  standard 
transformation  with  a q of  two  is  not  exact  for  turbulent 
flow;  but,  for  low  turbulence  levels  and  boundaries  with 
slight  curvature,  it  is  a good  approximation. 

The  difficulty  arises  when  transforming  the  term 
3t/3y.  Tetervin  shows  that  if  the  exponent  in  the  inte- 
grand of  the  transformation  on  x is  two,  then 

T Ri 

r = -r  • (48) 

T L 

This  corresponds  to  laminar  flow  or  flows  with  low  turbu- 
lence levels.  As  the  turbulence  increases  without  limit, 

1 - 1 , (49) 

T 

which  corresponds  to  an  exponent  of  unity  on  the  integrand 
of  the  x transformation.  In  this  instance,  it  is  a simple 
matter  to  determine  the  relationship  between  the  spreading 
parameters  also. 
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Consider  the  control  volume  of  Fig,  4,  page  34, 


y 2 

/ pu  Rdy 
yL 


then , 


y z, 

I pu  dy 


x ~ 7 o 

2 tt  f i Rdx  = 2uR  / pu^dy 


Now,  as  tq  is  independent  of  R for  high  turbulence  levels. 


R+R  o n 0 n ? 

T ( ) = R^PyU  -i-  / <Tdri  . 

0 2 . O u oA  n Pu 

L 


For  a two-dimensional  mixing  layer 


T = PuU2  i / f-  92dn 

J nL  U 


t 2R  - 
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"f  Rs  + Ro  aA 


Then 
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If  Rq  is  taken  at  the  inviscid  boundary  and  FL  is  con- 
sidered a good  approximation  to  R,  then 


2R. 

= ( - 1 ) o , 

Rs  + Ri 


C56) 


which  is  exactly  Eq.  (44).  The  Mangier  transformation  for 
this  case  gives 


or 


Ly 


f 

o 


X 


R . dx 
i 


(57) 


- a < 


2Ri 


Rs  + Ri 


) (J) 


(58) 


for  a conical  boundary. 

Thus,  the  modification  to  a given  by  Eq.  (44)  can 
be  interpreted  as  the  Mangier  transformation  for  highly 
turbulent  flow. 

For  slightly  turbulent  flow  and  an  exponent  of  two. 


n 


3R . R 
1 s 


Rt  + R.  R + R 

1 IS  5 


(59) 


where  a reference  length  of  Rg  was  used  to  be  consistent 
with  Eq.  (57).  That  is,  for  the  case  of  an  inviscid 
boundary' parallel  to  the  axis  of  symmetry  where  = Rg , 
there  is  no  distortion  of  the  spreading  parameter. 
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Dixon,  et  al . [30]  also  consider  to  be  appropriate  as 
the  reference  length. 

XI.  Inviscid  Field 


Basic  Method  of  Characteristics 

The  general  inviscid  flow  field  is  developed  using 
the  classical  method  of  characteristics.  A presentation  of 
the  details  of  the  method. is  given  by  Fox  [31] . An  inter- 
esting der Nation  of  the  method  is  given  in  Appendix  C. 
Starting  with  the  axially  symmetric,  inviscid  conservation 
equations  of  mass,  axial  and  radial  momentum,  and  energy, 
as 


pv 


3pv 

3r 


+ r 


i£B  = o 

3z  u 


and 


3u  + 5p  ' 
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(60) 
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3H 
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0 . 


With  the ' definition  of  total  enthalpy  as 


H = j ( u2  + v2  ) + h , (61) 

the  well-known  ordinary  differential  equations  of  the 
method  are  derived  as 

udr  - vdz  = 0 (62) 
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and 

( - ai  ) dr^  - 2uvdrdz  + { - a‘  ) dz^  = 0 (63) 

for  the  characteristic  surfaces.  The  compatibility 
relations  along  the  first  surface  which  is  a streamline  are 


and 


dp/p  + udu  + vdv  = 0 


(64) 


dh  + udu  + vdv  = 0 . 


(65) 


The  compatibility  relations  along  the  second  set  of  surfaces, 
which  are  the  left  and  right -running  characteristics,  are 

2 

(udr  - vdz) (udz  - vdu)  + (udr  - vdz)  v/r 

+ dp/p  (vdr  + udz)  = 0 (66) 

The  equations  are  put  in  finite-difference  form  and  inte- 
grated numerically  subject  to  the  conditions  at  the  plume 
boundary.  As  the  plume  boundary  is  a free  surface,  the  con- 
ditions are  : 

1.  The  surface  is  a streamline. 

2.  The  pressure  is  constant  along  the  surface. 

Also,  there  are  particular  features  of  the  flow  field 

that  require  special  treatment,  as  detailed  below. 

Embedded  Shock 

The  embedded  shock  within  the  plume  is  treated  in 
detail  in  [31];  consequently,  only  an  overview  will  be 
presented  here. 
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The  shock  wave  occurs  in  the  solution  field  when 
characteristic  curves  of  the  same  sense  coalesce  and  over- 
lap one  another.  Double- valued  points  occur.  A shock 
discontinuity  is  inserted  to  match  the  low-pressure  core 
region  to  the  high-pressure  boundary  region.  The  conser- 
vation equations  across  the  shock  reduce  to 


P2u0  sin  b + P-?v2  cos  ^ = P]_ui  sin  & + cos  ® = Cy  » 

CjU?  + p7  sin  b = sin  3 , 

Clv2  + p2  cos  S “ Clvl  + P1  sin  ® * 


(67) 


and 


, . U2  , V2  . U1  V1 
h2*V  T -S*  — *T 


where  the  subscript  "1"  refers  to  known  upstream  conditions, 
the  "2'’  refers  to  the  downstream  conditions  to  be  found, 
and  6 is  the  angle  between  the  shock  and  the  axial 
direction,  measured  clockwise  positive. 

Corner  Expansion 

At  the  corners  encountered  both  by  the  plume  and 
the  outer  stream,  the  flow  must  turn  abruptly  to  match  the 
base  pressure.  In  an  inviscid  flow,  the  flow  properties  at 
the  point  of  expansion  are  not  single  valued;  that  is,  the 
field  at  that  point  is  assumed  to  experience  a continuous 
change  in  pressure  from  the  pressure  existing  before  the 
turn  to  the  pressure  of  the  base  region.  The  required 
relations,  valid  at  a point,  are  called  the  Prandtl-Meyer 
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relations.  They  can  be  shown  to  be  the  limiting  forms  of 
the  compatibility  relations  along  both  a characteristic  and 
a streamline  as  the  coordinate  differentials,  dr  and  dz,  go 
to  zero.  They  are 


and 


du  _ 
dp 


2 2 
tr  + vL 


(68) 


dv 

dp 


- - /P^l 
_e p 


U‘ 


+ v^ 


(69) 


Thus,  the  velocities  can  be  determined  as  a function  of 
pressure  at  a point. 


Recompression  Shocks  at  the  Intersection  of  the  Streams 

Figure  6 is  a more  detailed  representation  of  the 
recompression  region  resulting  from  the  intersection  of  the 
two  streams.  The  downstream  side  of  both  shock  waves  must 
match  flow  angle  and  pressure  along  the  slipline  separating 
the  two  flows.  Therefore, 


and 
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These  two  relations  coupled  with  two  sets  of  shock  rela- 
tions, Eq . (67),  are  sufficient  to  solve  for  the  flow 
parameters  along  the  slipline  downstream  of  the  shock. 

Thus,  given  the  initial  conditions  upstream  of  the 
corners  for  both  streams  and  the  pressure  in  the  base 
region,  the  inviscid  boundaries'  are  developed  in  a marching 
scheme  using  the  method  of  characteristics  until  the  inter- 
section point  writh  the  recompression  region  is  determined. 


49 


A EDC-TR  -79-46 


CHAPTER  III 

SOLUTION  PROCEDURES 

I.  Introduction 

In  Chapter  II  the  overall  theory  for  the  sepa- 
ration, mixing,  and  compression  of  turbulent  viscous  layers 
was  presented  along  with  that  of  the  corresponding  inviscid 
field.  It  is  the  purpose  of  the  following  sections  to  show 
in  an  unambiguous  way  the  solution  procedure  used  to 
determine  the  base  properties  when  two  different  supersonic 
streams  (the  air  stream  and  the  rocket  exhaust  in  the  most 
general  case)  separate,  mix  with  the  base  gas,  and 
recompress  to  the  downstream  pressure.  The  principal 
assumption  employed  is  that  the  two  layers  do  not  interact 
except  as  they  both  impact  and  turn  upon  an  imaginary  slip- 
line. The  slipline  is  determined  completely  by  the 
inviscid  computations  using  the  method  of  characteristics. 
The  mixing  and  recompression  processes  of  the  two  mixing 
layers  are  assumed  to  occur  independently  of  one  another , 
having  only  the  base  pressure  in  common.  The  two  mixing 
fields  are  coupled  only  through  the  global  conservation 
relations.  All  of  the  terms  of  the  global  relations, 
except  the  bleed  terms,  are  evaluated  at  the  inviscid  in- 
tersection. The  initial  boundary  layer  is  assumed  to  ex- 
perience no  distortion  from  the  corner  expansion,  process. 
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II.  General  Approach 

The  primary  system  of  equations  to  be  solved  is 
the  conservation  of  mass 


'D  yD  yS  yS 

2tt  / puRdy  ] + 2tt  / puRdy  ] = m+  2tt  / puRdy  ] + 2ir  / puRdy  ] 


SI 
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S2 


SI  y. 


S2 


(72) 


and  the  conservation  of  energy 
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with  the  unknowns  being  pB  and  Hg . It  is  convenient  to  use 
the  mass  fractions 
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and  the  Crocco  relation  between  species  distribution  and 
the  local  total  enthalpy  as 

H = kHx  + (1  - k)  Hb  (76) 

for  the  mixing  layer  of  Stream  1,  and 

H = kH2  + (1  - k)  Hb  (77) 

for  that  of  Stream  2,  to  reduce  the  energy  equation  to 

H0  = Hj  mf1  + H2  mf 2 + (1  - mf 2 - mf2)  H£  . (78) 

For  ease  of  computation,  Eqs.  (72),  (74),  and  (75)  are 
taken  as-  primary  with  the  unknowns  being  pB,  mf^,  and  mf2. 
With  no  bleed  flow,  they  reduce  to  two  equations  and  two 
unknowns  pg  and  mf ^ . 

They  are  solved  using  the  classical  Newton- Raphson 
technique  with  numerical  derivatives.  (See  [32].) 

Although  the  solution  of  the  system  is  simple,  the 
evaluation  of  the  terms  is  tedious.  It  would  be  unpro- 
ductive to  reproduce  here  the  many  calculations  involved  as 
they  use  no  novel  numerical  techniques. 

III.  Computational  Scheme 

To  begin,  the  problem  must  be  defined  through  a set 
of  geometric,  thermodynamic,  and  fluid  dynamic  parameters. 
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1.  Radius  and  half-angle  of  the  nozzle. 

2.  Radius  of  the  afterbody. 

Thermodynamic 

1.  Relative  mole-atoms  of  carbon,  hydrogen, 
nitrogen,  oxygen,  and  argon  in  Stream  1 and 
Stream  2 and  in  the  bleed  gas . 

2.  Total  enthalpies  and  pressures  of  the  streams. 

Fluid  Mechanical 

1.  Exit  static  pressure  of  nozzle. 

2.  Boundary  - layer  thickness  at  exit  of  nozzle. 

3.  Static  pressure  of  free  stream. 

4.  Boundary-layer  thickness  at  separation  on 
afterbody. 

Then  to  begin  computation,  the  first  guess  of  base  pressure 
and  mass  fractions  of  the  streams  in  the  base  must  be 
provided . 

Inviscid  Flow 

As  the  total  conditions  for  each  stream  are  known, 
the  entropy  for  each  stream  can  be  calculated.  At  constant 
entropy,  the  streams  are  expanded  to  the  base  pressure 
where  the  gas  properties  of  the  two  streams  are  determined. 
The  starting  lines  are  then  determined  for  the  method  of 
characteristics.  The  method  of  characteristics,  using  the 
frozen  equilibrium  gas  properties,  provides  the  entire 
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inviscid  flow  field.  Most  importantly,  it  allows  the 
length  of  the  inviscid  boundary  to  the  point  of  inter- 
section of  the  rocket  stream  and  the  free  stream  to  be 
determined.  This  information  is  required  for  entering  the 
viscous  flow  calculations. 

Viscous  Flow 

All  calculations  of  the  viscous  flow  are  made  at 
one  point  for  both  streams:  the  intersection  of  the 
inviscid  boundaries.  This  is  because  the  intersection  is 
on  the  slipline  dividing  the  two  streams ; the  slip  surface 
is  assumed  to  be  the  reaction  plane  upon  which  the  merging 
streams  impinge  and  divide  into  the  parts  which  are  turned 
back  into  the  base  and  those  which  are  recompressed  to  the 
downstream  pressure.  Thus,  the  properties  of  the  mixing 
layer  positioned  at  that  point  are  those  used  for  the 
evaluation  of  the  terms  in  the  primary  equations. 

The  spreading  parameter,  a,  is  calculated  from 
Eq.  (42).  As  this  spreading  parameter  is  based  on  fully 
developed  flow,  the  elemental  species  distribution  is 
assumed  to  follow  the  velocity  profile  from  the  Crocco 
relation.  That  is, 

k = u/U  . (79) 

(A  full  discussion  of  the  evaluation  of  a typical  integral 
is  provided  in  the  following  section.)  The  axisymmetric 
correction  to  a is  applied  using  either  Eq.  (58)  or 
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Eq.  (59),  depending  on  the  assumption  of  the  turbulence 
level.  Equations  (17)  and  (18)  can  now  be  solved  simul- 
taneously to  obtain  y and  C The  dividing  streamline 
location,  y^,  can  then  be  determined  from  Eq.  (1).  The 
calculations  are  repeated  for  each  mixing  layer. 

The  equations  of  the  recompression  analysis  must 
now  be  solved  iteratively  for  each  mixing  layer  to 
determine  the  location  of  the  stagnating  streamlines,  y^ . 
Thus,  Eqs.  (30),  (31),  (32),  (33),  and  (34)  must  be  solved 
together  with  geometrical  relations^  for  and  Xg  to 
determine  >'g  for  each  layer. 

All  the  parameters  necessary  for  evaluating  the 
primary  equations  are  now  available.  The  terms  are  evalu- 
ated and  the  residuals  determined.  The  Newton- Raphs on 
method  serves  to  predict  new  guesses  for  p^,  mf ^ , and  mf,, 
and  the  global  iteration  proceeds  to  convergence. 

Evaluation  of  Typical  Integral 
Consider  the  integral 

>'S 

I,  = f kpuRdy  . (80) 

yL 

The  local  mass  fraction  of  the  high-speed  stream  in  the 
local  mixture  at  a point  in  the  mixing  layer  is  defined  by 


■'‘These  relations  are  as  presented  in  [2]  and  are 
easily  derived. 
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Eq . (16).  Repeated  here,  it  is 


k = | [ 1 + erf  (Cmf  ] . (81) 

The  local  total  enthalpy  is 

H = kHu  + (1  - k)  Hb  . (82) 

The  velocity  at  any  given  y is  known  from  the  velocity 
profile,  Eq.  (13).  Also,  the  relative  moles  of  each 
element  are  determined  as 


Nj  = kNjcc  + (!  - k)  NjB  j=1>m  * CS3) 

Since  the  pressure  is  known  and  the  static  enthalpy  is 

h = H - u2/2  (84) 

the  mole  fractions  of  the  molecular  constituents  and  the 
temperature  are  then  determined  from  the  equilibrium 
chemistry;  then 


l C-w. 
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) . 


The  radius  to  any  point  in  the  mixing  zone  is 


(85) 


R - Rj_  + (y  - ym  ) cos  $ . (86) 

Thus,  all  elements  of  the  integrand  are  readily  found,  and 
the  integration  is  performed  using  a conventional  quad- 
rature. 
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CHAPTER  IV 

RESULTS  AND  DISCUSSION 
I . Introduction 

Checking  the  theory  with  experimental  data  began 
with  the  classical  two-dimensional,  backward- facing  step 
problem.  Quite  a large  amount  of  data  exists  for  this  con- 
figuration [8].  It  served  not  only  as  a check  on  the 
recorapress ion  model,  but  it  provided  critical  data  on  the 
effect  of  the  initial  boundary  layer. 

Davis  [33],  for  his  dissertation,  performed  some 
interesting  experiments  by  bleeding  hydrogen  gas  into  the 
base  region  of  the  backward-facing  step  and  taking  measure- 
ments both  with  and  without  burning.  This  provided  the 
information  to  check  the  effects  of  both  base  bleed  and 
chemistry. 

There  remained  only  to  find  data  for  the  configu- 
ration of  interest;  that  is,  an  axisymmetric  jet  embedded 
in  a blunt-based  afterbody.  The  Reid  and  Hastings  experi- 
ments [34]  were  ideal  for  this.  Their  experiments  involved 
exhausting  an  embedded  supersonic  jet  into  a supersonic 
stream  at  various  j et- to- stream  total  pressure  ratios. 

These  experiments  served  to  check  the  axisymmetric  approxi- 
mation. And  finally,  some  unpublished  data  from  an  AEDC 
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test  [35],  although  incomplete,  provided  a comparison  with 
a hot  rocket  exhausting  into  a stream  of  supersonic  air. 

II.  Backward- Facing  Step 


Isoenergetic  Flows 

The  data  as  compiled  by  [4]  are  presented  in 
Figs.  7,  8,  and  9 along  with  the  corresponding  predictions 
of  the  present  theory.  It  is  assumed  that  the  inviscid 
stream  approaching  the  backstep  is  uniform  and  extends 
without  interference  to  infinity  in  the  direction  trans- 
verse to  the  floxir. 

The  axially  symmetric  approximation  was  programmed; 
so  the  planar  problem  was  simulated  by  using  the  embedded 
nozzle  in  a blunt  - afterbody  configuration.  The  nozzle  was 
extended  in  radius,  in  comparison  to  the  height  of  the 
base,  until  the  axisymmetric  effect  was  negligible;  that 
is,  the  ratio  of  the  nozzle  radius  to  the  base  height  was 
50:1.  The  nozzle  half-angle  was  taken  as  zero,  and  the 
conditions  on  the  exit  plane  made  identical  to  those  of  the 
outer  flow. 

To  test  the  effect  of  the  momentum  thickness 
Reynolds  number,  the  correlation  given  graphically  by  [4] 
was  used.  It  is  reproduced  in  Fig.  TO. 

The  relatively  flat  trajectory  of  the  theoretical 
prediction  of  base  pressure  with  respect  to  momentum 
thickness  indicated  that  the  recompression  model  was 
behaving  correctly.  The  recompression  model  of  [2]  was 
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Figure  8.  Base  pressure  vs  momentum  thickness  of  initial 
boundary  layer  . (Moo  = 2.0]. 


A E □ C-TR  -79-4  6 


Figure  9.  Base  pressure  vs  momentum  thickness  of  initial 
boundary  layer  = 3.0). 
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based  on  an  inviscid  assumption.  It  is  shown  for  com- 
parison in  Fig.  8.  Data  scatter  is  typical  of  these 
experiments.  Part  of  this  is  attributed  to  wall  effects, 
and  part  is  due  to  the  different  histories  of  the  flows. 

(See  [15].) 

In  Fig.  11,  the  recompression  factor  of  Nash  [10] 
is  predicted  as  a function  of  the  boundary-layer  momentum 
thickness.  These  recompress  ion- factor  curves  are  those 
which  correspond  with  the  base  pressures  of  Figs.  7,  S,  and 
9.  As  Nash  did  not  take  into  account  the  initial  boundary 
layer,  comparisons  are  difficult.  However,  he  stated  that  a 
recompression  factor  of  0.35  correlated  data  well  for 
supersonic  Mach  numbers.  In  the  present  case,  say  at  a 
Mach  number  of  2.0,  a recompression  factor  of  0,35  corre- 
sponds to  a 9/h  of  about  0.038,  which  is  quite  typical. 

Perhaps  the  most  severe  test  of  the  present  theory 
was  afforded  by  the  data  [36]  presented  in  Fig.  12.  Here, 
computations  were  made  to  boundary- layer  thicknesses  which 
were  approximately  nine  times  the  step  height.  After 
approximately  a 9/h  of  0.8,  convergence  became  very  slow; 
and  finally,  convergence  could  not  be  achieved  at  all. 

At  a Mach  number  of  1.5,  Fig.  7,  the  predictions 
fall  low  in  the  data.  However,  Fig.  13  compares  the  pre- 
dictions with  data  from  Tanner  [9].  Although  he  gives  no 
initial  boundary- layer  information,  it  was  found  that  a 
boundary- layer  momentum  thickness  of  0.04  agreed  almost 
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Figure  12.  Base-pressure  predictions  at  very  large  initial 
boundary- layer  momentum  thicknesses. 
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Figure  13,  Bas e -pressure  coefficient  as  a function  of 
Mach  number  (Tanner's  data). 
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exactly  with  his  theory,  which  he  claims  to  be  a good  pre- 
diction of  base  pressures  in  the  low  supersonic  range. 

From  these  results,  it  was  concluded  that  the 
theory  worked  quite  well  for  cold,  planar  flows  writhout  any 
base  bleed.  The  next  test  would  be  the  predicting  of  base 
pressures  with  base  bleed  and  with  burning. 

Base  Bleed  and  Burning 

Davis  [33]  performed  the  following  experiments.  A 
blunt-based  afterbody  was  immersed  in  an  air  stream  flowing 
at  a Mach  number  of  2.0.  Hydrogen  was  bled  into  the  base 
region  at  various  flow  rates.  The  base-pressure  change  was 
plotted  against  the  flow  rates  producing  the  lower  data 
curve  of  Fig.  14.  Then,  a pilot  flame  was  ignited  and  the 
bleed  flow  was  adjusted  to  the  desired  rate.  The  flame  was 
extinguished  when  the  wake  burning  was  self-sustaining. 

The  base  temperature  and  pressure  were  plotted  producing 
the  upper  curve  of  Fig,  14  and  the  data  curve  of  Fig.  15. 

Computations  were  made  to  simulate  the  experiments. 
The  results  are  plotted  with  the  data  curves  of  Figs.  14 
and  15.  As  Davis  obtained  no  initial  boundary-layer  data, 
the  boundary- layer  thickness  of  the  theory  was  adjusted 
until  the  computed  base  pressure  with  no  bleed  matched  the 
zero  bleed  flow  point  of  the  data.  With  a oq  of  12.0,  this 
pressure  was  obtained  at  a 0/h  value  of  0.019.  At  a oo  of 
9.0,  a 0/h  value  of  0.04  was  required.  When  burning  was 
initiated  in  the  computations  at  a oQ  of  9.0,  no  solutions 
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Figure  14.  Effect  on  the  base  pressure  of  a hydrogen  bleed 
both  with  and  without  burning. 
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were  found,  The  curves  presented  are  those  calculated  at  a 

oq  of  12.0,  It  is  not  clear  that  any  physical  significance 

can  be  attached  to  conditions  with  a 0 of  9.0,  for  even 

o 

with  the  computations  at  a aQ  of  12.0  after  a bleed  rate  ■ 
parameter,  H,  of  6.0  x 10  ^ , convergence  became  very  slow; 
and  then  finally,  no  convergence  could  be  obtained.  Per- 
haps all  that  can  be  reasonably  said  is  that  the  limits  of 
the  theory  are  being  observed  and  nothing  more. 

The  pressure  prediction  is  high,  but  acceptable. 

The  temperature  prediction  is  also  high;  however,  the  loca- 
tion of  the  peak  is  well  predicted,  which  indicates  that 
the  basic  flow  mechanism  is  well  modeled.  It  is  inter- 
esting to  observe  that  in  Davis'  theoretical  approach,  he 
input  the  base  pressure  as  knowns  to  obtain  the  base 
temperatures.  As  pointed  out  in  Chapter  II,  Section  VII, 
his  flame-sheet  model  disallows  oxidizer  in  the  base  region, 
which  accounts  for  his  lower  temperature  predictions,  for 
the  peak  temperatures  in  the  mixing  layer  of  approximately 
4300°R  are  predicted  by  his  method,  which  matches  the  peak 
temperatures  of  the  present  theory.  That  is,  if  oxidizer 
is  in  the  base,  at  some  bleed  rate  stoichiometric  con- 
ditions must  occur.  It  is  perhaps  also  worthy  of  note  that 
the  adiabatic  equilibrium  conditions  as  assumed  by  the 
theory  are  rarely  achievable  in  real  experiments  of  this 
type  due  to  the  effects  of  wall  cooling,  thermocouple 
errors,  and  combustion  inefficiency. 
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The  curves  of  Fig.  16  are  presented  to  show 
temperature  profiles  as  predicted  by  the  theory.  No  data 
are  presented  as  Davis’  recordings  across  the  mixing  layer 
were  too  sparse  to  encompass  this  area. 

In  conclusion,  here,  it  must  be  said  that  the  pre- 
dicted shapes  of  the  curves,  particularly  the  base  pressure 
curves,  show  that  the  basic  physics  of  this  configuration 
is  well  predicted;  the  theory  requires  only  a more 
realistic  chemistry  model. 

III.  Correction  for  Axially  Symmetric  Effects 

As  previously  pointed  out,  the  theory  is  two- 
dimensional.  However,  the  corrections  for  axisymmetric 
effects  were  employed  as  detailed  in  Chapter  II,  Section  X. 
Reid  and  Hastings  [34]  recorded  the  effect  of  total 
pressure  ratio  as  between  the  embedded  jet  and  the  external 
stream,  both  of  which  were  at  a Mach  number  of  2.0.  Com- 
plete initial  boundary- layer  information  was  provided  which 
allowed  the  best  check  yet  of  the  present  theoretical 
approach.  For  the  configuration  where  the  outside  radius 
is  of  a ratio  of  1.67  to  the  nozzle  radius,  the  following 
information  was  provided: 
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Figure  16.  (Continued) 
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Nozzle 


External  Body 


1.17  in. , 

0.0037  in . , 

2,78  in  the  base  plane. 

1.95  in . , 

0 . 0165  in.  , 

2.97  in  the  base  plane. 


A power  law  exponent  of  1/9  and  boundary- layer 
thicknesses,  = 0,048  and  6,/Rj  = 0.212,  gave  the 

correct  momentum  thicknesses.  Boundary- layer  information 
was  not  given  on  the  other  nozzles;  but  since  the  aspect 
ratios  of  the  nozzles  were  approximately  equal,  the  value 
of  0.048  for  S^/R^  was  used  for  the  others  as  well.  The 
external  momentum  thicknesses  were  the  same  for  all  con- 
figurations . 

Using  an  incompressible  reference  spreading 
parameter,  a , of  9.0,  and  a as  derived  in  Chapter  II, 
Section  X,  the  results  presented  in  Fig.  17  were  obtained. 

The  prediction  at  R?/R^  equal  to  1.25  was  particu- 
larly welcomed  as  it  further  confirmed  the  basic  theory  as 
this  case  is  close  to  being  two-dimensional.  As  the  radius 
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and  Hastings  [34] 
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Figure  17.  Relationship  between  base  pressure  and  the 
jet-to-free-stream  total  pressure  ratio  of 
axially  symmetric  jet. 
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ratio  is  increased  and  the  axisymmetric  effects  become  more 
apparent,  the  approximation  is  not  quite  as  good,  but  still 
acceptable.  For  these  cases,  the  conventional  Mangier 
transformation  which  produces  the  spreading  parameter  of 
Eq.  (59)  gives  the  best  results;  however,  the  of  Eq.  (59) 
appears  to  better  follow  the  shape  of  the  data  curve. 

The  last  axisymmetric  configuration  attempted  was 
at  once  the  most  complex  and  the  most  incomplete.  Its 
complexity  stems  from  the  fact  that  it  is  a hot-rocket 
exhaust  (y  = 1.26)  embedded  in  a relatively  large  base  sub- 
merged in  a supersonic  air  stream.  It  is  incomplete 
because  the  initial  boundary  layer  is  defined  only  by  its 
thickness,  and  the  base  pressure  is  the  only  parameter  that 
was  measured.  Nevertheless,  a series  of  runs  was  made  at 
external  stream  Mach  numbers  of  2.05  and  3.00,  the  results 
of  which  are  shown  in  Fig.  18. 

The  external  boundary- layer  thickness  was  main- 
tained at  67  percent  of  the  nozzle  radius,  and  the 
different  momentum  thicknesses  resulted  from  varying  the 
shape  of  the  profile.  The  results  at  the  lower  Mach  number 
are  reasonable  although  a relatively  large  momentum 
thickness  is  required  to  match  the  recorded  base  pressure. 

In  all  cases  the  boundary- layer  momentum  thickness  of  the 
jet  was  maintained  at  approximately  2 percent  of  the  base 
thickness . 
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Figure  18. 


Calculated  base  temperatures  and  pressures  for 
a hot  rocket  at  two  external  stream  Mach 
numbers  with  varying  external  boundary- layer 
momentum  thicknesses. 
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The  results  at  the  higher  Mach  number  were  dis- 
appointing in  that  near  the  recorded  base  pressure 
convergence  was  not  achieved.  What  is  apparent,  however, 
is  the  effect  of  the  external  boundary  layer  on  the  temper- 
ature and  pressure  of  the  base. 

In  both  Mach  number  cases,  the  temperature  drops 
writh  increasing  external  boundary- layer  momentum  thickness. 
This  is  because  more  of  the  external  stream  and  less  of  the 
rocket  exhaust  are  captured  in  the  base  region.  Because  of 
the  much  larger  thickness  of  the  outer  boundary  layer,  how- 
ever, it  tends  to  be  the  overriding  mechanism  affecting  the 
base  pressure. 

All  calculations  for  this  case  were  made  assuming 
frozen  chemistry  as  equilibrium  chemistry  results  were 
quite  unrealistic.  They  indicated  that  burning  was 
occurring,  which  served  to  drive  pressures  and  temperatures 
much  higher  than  experiment  would  indicate. 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

A generalized  supersonic  base-flow  theory  has  been 
developed  for  the  near-wake  region  behind  bluff-based,  two- 
dimensional  and  axially  symmetric  bodies  in  the  presence  of 
a coaxial,  supersonic  jet.  The  theory  works  well. 

The  correlations  presented  show  the  theory  to  be 
viable  and  physically  realistic.  In  particular,  a new 
recompression  model  was  able  to  generalize  the  mixing 
analysis  and  extend  the  overall  theory  to  encompass  many 
complex  and  diverse  flow  configurations,  including  those 
with  very  thick  boundary  layers . The  new  species  distri- 
bution model  allowed  flows  with  hydrogen  bleed  to  be 
predicted  as  well  as  flows  with  large  total  temperature 
differences  between  the  base  and  the  high-speed  streams. 

As  the  mechanism  v/as  available,  a chemical  equilibrium 
analysis  was  coupled  to  the  theory.  It  gave  an  excellent 
prediction  of  the  point  at  which  the  base  bleed  rate  would 
produce  the  maximum  base  temperature.  It  also  gave  a good 
prediction  of  the  effect  of  burning  on  the  base  pressure. 

The  Mangier  transformation  was  used  to  develop  an 
axially  symmetric  correction  to  the  two-dimensional 
spreading  parameter  for  highly  turbulent  mixing.  The 
transformation  also  was  used  to  provide  an  axi symmetric 
correction  to  the  spreading  parameter  for  low -turbulence 
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and  laminar  flows.  However,  interesting  and  useful  as 
these  corrections  were,  they  further  amplified  an  under- 
lying problem:  the  lack  of  a good  model  for  the  spreading 
parameter . 

Indeed,  perhaps  the  most  important  result  of  this 
investigation  was  the  realization  that  the  modeling  of  the 
spreading  parameter  needs  further  development.  There  is  a 
tendency  in  the  literature  to  treat  the  spreading  parameter 
as  an  off-the-shelf  item,  important  only  to  the  extent  that 
some  numerical  value  must  be  chosen.  This  is  understand- 
able, of  course,  when  the  heavy  hand  of  empiricism  must  be 
used  to  nudge  an  unsatisfactory  solution  into  the  area  of 
reality,  as  is  necessary  when  a recompression  factor  is 
used,  and  the  initial  boundary  layer  is  ignored.  However, 
as  the  basic  theory  becomes  more  realistic  and  inclusive, 
the  necessity  of  a more  realistic  spreading  parameter 
becomes  apparent.  Specifically,  the  spreading  parameter 
must  mimic  reality  to  the  extent  that  it  develops  as  the 
mixing  develops  and  is  asymptotic  to  the  fully  developed 
value,  taking  into  account  the  history  of  the  separating 
boundary  layers. 

Although  turbulent  flows  were  of  primary  interest 
in  the  present  analysis,  appropriate  modeling  of  the 
spreading  parameter  would  allow  the  theory  to  be  extended 
to  laminar  mixing.  This  becomes  important  in  low-density 
flows  . 
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Perhaps  the  weakest  point  in  the  present  theory  is 
the  lack  of  a good  model  for  axially  symmetric  mixing.  The 
assumption  of  planar  mixing,  while  giving  acceptable  pre- 
dictions for  the  most  part,  did  begin  to  fail  at  afterbody - 
to-nozzle  radius  ratios  greater  than  1.67. 

Again,  in  overview,  it  must  be  said  that  the 
present  theory  works  well.  It  predicted  the  effect  of  the 
initial  boundary  layer  in  planar  flow  as  well  or  better 
than  more  sophisticated  analyses  up  to  very  high  momentum 
thicknesses.  The  theory  predicted  well  the  effect  of 
hydrogen  bleed  both  with  and  without  burning.  Predictions 
of  axially  symmetric  flows  with  coaxial  cold  jets  were 
excellent  for  small  bases  and  quite  acceptable  for  the 
larger  bases.  It  is  concluded  that  the  present  theory  is 
at  least  the  equal  of  the  Crocco-Lees  approach  to  simple 
flows,  and  it  is  superior  when  applied  to  complex  flows 
such  as  the  two-mixing-layer  problem. 

The  extension  to  flows  with  nonequilibrium 
chemistry  should  be  a fertile  ground  for  further  work.  By 
coupling  a rate  equation  to  the  system  and  adding  the  non- 
equilibrium term  to  the  thermodynamic  property  relations, 
the  present  theory  should  be  able  to  predict  the  effect  of 
nonequilibrium  on  the  near  wake. 
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APPENDIX  A 

DETERMINATION  OF  EQUILIBRIUM  GAS  COMPOSITION 
BY  SUCCESSIVE  APPROXIMATIONS 


It  is  the  purpose  here  to  show  by  example  the 
method  of  Weinberg  [37]  as  developed  by  Osgerby  and 
Rhodes  [38]  to  determine  the  equilibrium  gas  composition 
(mole  fractions)  of  a mixture  as  a function  of  pressure, 
temperature,  and  the  relative  number  of  the  various 
elements  presents 

The  C,H,Q  case  is  used  to  illustrate  the  method. 
The  equilibria  are 


co2  : CO  * \ 02  , 

H2°  H2  + 2 °2  ’ 
h2o  ; [ h2  * oh  , 


1 

1 


) 


i 


(A.l) 


and  the  atom  conservation  equations  are 


p0  “ 2pC02  + pCO  + Ph20  + 2p02  + P0H  + P0  ' 
PH  = 2pH20  + 2pH2  + P0H  + PH  ’ 

PC  = PC02  + PC0  » 
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s 


/ 

where  P is  the  fictitious  pressure  of  the  particular 
element . 

The  conservation  equations  are  rewritten  to  elimi- 
nate the  fictitious  pressures  by  using  the  fact  that 


3 N 
H J H 


HrC  , 


(A. 3) 


0 _ 0 


0rC  , 


(A. 4) 


and 


pco2  + Pco  + Ph2o  + Po2+Poh  + Po  + ph2+Ph 


= p 


(A.  5) 


The  new  conservation  relations  are,  then, 


F"2pCo/PcO  + PH20+2Po2+P0H  + Po‘0rC(PC02+PCO)=0  ’ CA-6] 


G = 2pH20  + 2pH2  + P0H  + PH  " HC  ^ PC02  + PC0  ^ = 0 ' 


(A.  7) 


and 


L-PCO,*PCO*%,0*PO  *PcH^0*>'H,+PH-p'0  • (A’8) 


From  Eq.  (A.l),  the  equilibrium  relations  may  be 
written  as 


>C0  = K! 


Xo, 


yPr 


H, 


K. 


’h2o=  k2  ^ 


P0H  " K-  P'0 


2 2 
(A.  9) 


PH  = k4  ^ and  p0_K5  f 
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where  K is  the  equilibrium  constant  defined  as 


K = exp  [ 


(6j  -«j) 


RT 


(A. 10) 


where  the  chemical  potentials,  y°  (or  Gibbs  free  energy 
per  mole)  , are  as  given  by  [20]  . 

Since  Eq . (A. 9)  can  be  arranged  such  that  Pq,  Pq^» 
Pco,  pH,  and  pR  Q are  explicit  functions  of  pH2>  Pq2 , and 
PgO^'  they  may  be  substituted  into  Eqs . (A, 6)  to  (A. 8). 
This  gives  three  equations  for  the  three  unknowns,  p^, 

P02’  and  PC02- 

A Newton - Raphson  procedure  is  then  applied  which 
converges  very  rapidly.  That  is,  Eqs.  (A. 6)  to  (A. 8)  are 
linearized  [32]  as 


6F 

<5G 

5L 


^7  SPh2  + 5Po7  6p°2  * 4o,  SPc02  ’ 


8t  6pu  * «P 


”-hz  ■ SPoT  -o2  3Pco  -co2 


8L  p 3L  r ^ 3L 

5Pu  6P^  + ~ 6P 


wfV  Wo,  “rcoz 


(A. 11) 


A first  guess  of  pH^,  pQ  , and  pCQ  gives  the  initial 
residuals,  F,  G,  and  L,  and  also  the  values  of  the  partial 
derivatives.  These  three  equations  (Eq.  (A. 11))  are 
linear  equations  for  6pH  6pg  , and  Spgg  with  the  new  F, 
G,  and  L being  set  to  zero;  that  is, 
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SF  = -F  , 

5G  = -G  , (A. 12) 

5L  = -L  . 


The  new  guesses  of  the  unknowns  are  then 


^new  ^initial 


= 6p 


(A. 13) 


The  new  partial  pressures  become  the  first  guesses,  and  the 
procedure  is  continued  until  the  residuals  are  satis- 
factorily small. 
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APPENDIX  B 

A SUGGESTED  RESOLUTION  OF  THE  PROBLEM  OF  THE  PROFILE 
SHIFT  WHEN  CONSIDERING  AXIALLY  SYMMETRIC  MIXING 

In  the  following,  a method  is  suggested  for 
determining  y and  for  axially  symmetric  mixing.  That 

calculations  based  on  this  approach  are  not  presented  is 
because  the  approach  was  not  apparent  until  late  in  the 
life  of  the  project.  However,  it  will  be  used  as  a basis 
for  further  refinement  of  the  theory  to  be  taken  up  at  a 
later  date. 

Equations  (17)  and  (18)  of  the  two-dimensional 
theory  are  repeated  here  as 

yu  yU+ym 

/ p u dy  ] = / kpudy  ] (B . 1) 

o 5 S x=o  yL  x=L 

and 

XU  yU+ym 

/ p u 2 dy  ] - / PU2dy  ] CB- 2) 

o S 5 x=o  yL  X=SL 

with 

k = i [ l + erf  ( Cm£n  ) 1 • tB-3^ 

For  two-dimensional  mixing,  Eqs . (B.l)  and  (B.2)  are 

sufficient  for  determining  both  Cmf  and  ym. 
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In  words,  these  equations  say  that  when  the  mixing 
profile  is  localized  such  that  y is  the  distance  from  the 
origin  to  the  inviscid  boundary,  at  x = £,  then  the  mass 
and  momentum  of  the  Initial  profile  out  to  the  streamline, 
y^,  are  preserved,  The  streamline,  y^,  is  determined  at 
x = £ such  that  9 = 0.99998895,  In  two-dimensional  flow, 
the  streamline,  y^ , can  be  traced,  parallel  to  the  inviscid 
boundary,  to  the  beginning  of  mixing.  With  axially 
symmetric  mixing,  y^  is  not  parallel  to  the  inviscid 
boundary.  At  this  point,  for  axially  symmetric  mixing, 
there  is  not  enough  information  to  solve  the  problem.  . 

■It  is  important,  also,  at  this  point  to  stress 
certain  elements  of  the  theory. 

1.  The  boundary  layer  is  considered  as  part  of  the 
inviscid  flow;  that  is,  the  boundary  layer  is 
an  inviscid,  rotational  profile  and  no  more. 

2.  Equations  (B.l)  and  (B.2)  are  actually  rela- 
tions between  the  inviscid  flow,  the  left-hand 
side,  and  the  viscous  mixing,  the  right-hand 

5 ide  , 

With  this  in  mind,  Eqs.  (B.l)  and  (B.2)  for  two- 
dimensional  mixing  may  also  be  formulated  as 


ym+yU  yU+ym 

/ p u dy  ] = / kpudy  ] (B.-4) 

5 5 z-1  yL  X-I 
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and 

ym+yU  yU+ym 

/ P„u  2dy  ] = / pu2dy  ] . (B.5) 

ym  X=*  yL  X=* 

That  these  relations  are  equivalent  to  Eqs . (B.l) 
and  (B.2)  requires  only  the  observation  that  an  inviscid 
profile  in  two-dimensional  flow,  subject  to  no  pressure 
gradients  or  dissipative  mechanisms,  remains  intact.  In 
other  words,  Eqs,  (B.4)  and  (B.5)  say  that,  locally,  the. 
mass  and  momentum  rates  of  the  viscous  layer  are  equal  to 
the  mass  and  momentum  rates  of  the  inviscid  field  from  its 
boundary  to  the  same  y location,  at  the  same  x station. 

With  these  observations  in  mind,  the  axially 
symmetric  equivalent  of  Eqs.  (B.4)  and  (B.5)  may  be  formu- 
lated as 

ym+yU  yU+ym 

Ztt  / (pu)  *Rdy  ) = 2tt  / kpuRdy  ] (B.6) 

y x=2  y, 

and 

VyU  yU+ym 

2u  / (pu2)  *Rdy  ] = 2tt  f pu2Rdy  ] . (B.7) 

y x=Jt  y, 

The  asterisk  used  above  and  in  the  following  serves 
to  keep  track  of  the  fact  that  in  the  axially  symmetric 
case,  the  boundary  layer  thickness  and  shape  are  not  the 
same  at  x = 1 as  they  are  at  x = o,  That  is,  now  the 
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problem  is  complicated  by  the  fact  that  if 

1 

us/U  = ( y/6  )n  (B  . 8) 


is  the  initial  velocity  profile,  at  x = II 

1 

u/U  = (y/6*  )n  , (B . 9) 

where 

(B.10) 

and 

n*  + a . (B . 11) 

The  parameters,  6*  and  n“,  may  be  found  as  follows, 
Consider  that  the  following  relations  hold  for  the  inviscid 
flow  at  constant  pressure  and  with  no  dissipative 
mechanisms , 


_ . y +5*  , , * 

6 1/n  'm  1/n 

2tt  / p U ( y/6  ) . Rdy  ] = 2ir  / pU  (y/6*)  Rdy  ] 

o s x=o  y x=Z 

'm 

and 

6 2/n  Vf  , 2/n* 

2tt  / p U2  ( y/5  ) Rdy  ] = 2ir  / pU2  (y/6  ) Rdy]  . 
o s x=o  y_  x=l 


(B . 12) 


(B  . 1 3) 


Thus,  there  now  exist  sufficient  relations  to  solve 
the  axially  symmetric  problem.  The  four  relations,  then, 
Eqs . (B.6),  (B . 7)  , (B.12),  and  (B.13)  along  with  Eq . (B.3) 
serve  to  determine  the  four  unknowns , 6*,  n*  , Cmf,  and  ym- 
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It  is  realized,  of  course,  that  as  these  compu- 
tations will  be  embedded  in  the  overall  base -pressure 
iteration,  they  further  degrade  the  efficiency  of  the 
global  calculations. 
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APPENDIX  C 

DEVELOPMENT  OF  THE  CHARACTERISTIC  AND 
COMPATIBILITY  EQUATIONS 

Consider  the  conservation  equations  for  mass,  2 and 
r,  momentum  and  energy  as 


with 


and 


pv  + r + r = 0 , 


dr 


pu 


3u  ^ 3u  . 3p  n 
- — + pv  — + = 0 


3 2 


9r  3z 


3v  3v  9p  n 

pu  II  + pv  w + w 3 0 * 


3H  x 3H  n 

pu  _ + pV  __  = 0 , 


H = i ( uZ  + v2  ) + h 


p = f ( p , h ) . 


(C.l) 

(C.2) 

(C.3) 

CC.4) 

CC.5) 

(C.6) 


These  constitute  the  set  of  equations  for  u,  v,  p,  and  h. 
Now,  using 

— ft  i\  . M. 

(C.7) 


3 P _ 3p_  £_P.  + juo  _3_h 
r 


az  •*  3 p 3 z ah  3 z 


and 


3p  ' 3p  _3p_  + 3_P  JL!i 

Tr  \ 3p  3r  &h  3r  ’ 


(C.8) 
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Eqs • (C. 


where 


1)  through  (C.4)  may  be  cast  in  the  form 


(C . 9) 


pu 


n ii  3p  u 3p 

0 u -5p  u aH 


C pu 


pu  puv 


pu 


(C.10) 


B = 


pv 


P 

0 

0 pv 
puv  pv2 


v 


_3p 

3p 

■0 


3p 

3h 

0 

0 

pv 


Cc.ii) 


c = 


2l 

r 


0 

0 

0 


(C . 12) 
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and 


u 

V 

p 

h 


(C.13) 


It  is  desired  to  find  curves  in  the  region  of 
interest  along  which  the  set  of  partial  differential 
equations  reduce  to  ordinary  differential  equations.  Con- 
sider the  definition  of  a total  differential 


3W  a,  * 9W 

dW  = -RT-  dz  + -^r  dr  . 


3z 


9r 


(C.14) 


Along  a given  curve,  dr  and  dz  are  no  longer  independent, 
so  that 


dW  = 9W  + 9W  dr 
Sz  9 z 9r  3z 


8W 


(C.15) 


Using  Eq.  (C.15)  to  eliminate  ^ , 

dz  9r  cTz  ^ “ 9r 


A(f  - *E-o 


CC.16) 


is  obtained. 

Taking  the  remaining  partial  derivative  to  the 
right-hand  side  gives 


a|+  C -(gA-B)|| 


CC.17) 
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Now,  it  is  observed  that  if  a vector  were  found  that  was 
orthogonal  to  the  vectors  of  Eq . (C.17)  and  the  inner 
product  were  taken,  Eq.  (C.17)  would  appear  as 

( A g|  + C ) • A = 0 , (C.18) 

which  is  the  form  that  is  sought;  that  is,  Eq.  (C.18)  is  an 
ordinary  differential  equation.  Thus,  associated  with  each 
a vector  is  found  such  that 

dz  ’ 

( A 4^  - B ) • A = 0 . (C.19) 

K dz  or 

To  solve  for  such  a vector,  a property  of  the  scalar 
product  is  used  to  put  Eq . (C.19)  in  the  form 

M . (dr  . B)TA=0  , (C.20) 

3r  dz 

where  the  superscript,  T,  means  transpose.  Now,  since  ^ 
is  assumed  to  be  nonzero,  it  is  sufficient  that 

(A  - b/  A - 0 , (C.21) 


where  0 is  the  null  vector. 

From  linear  algebra,  for  there  to  exist  nontrivial 

solutions  to  a homogeneous  set  of  equations,  the  determi- 
nant of  the  coefficient  matrix  must  be  zero.  Thus, 


Det 


(C . 22) 
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This  condition  forces  specific  conditions  on  af  as  u is  a 

A Y* 

polynomial  in  -ji . These  are  called  the  characteristic 
equations.  For  this  case,  they  are 


c&>  •(&) 


31 


u 


(C . 23) 


and 


uv  + a 


47 


+ v 


2 . 


CC.24) 


, dr,  _ uv 
=“ 


i/u^  + 

2 2 
- a*1 


(C.25) 


dr 

These  relations  are  substituted  one  at  a time  back  into 

the  homogeneous  set  (Eq . (C.21)),  and  the  components  of 

four  different  vectors  are  found  corresponding  to  each  Sjy. 

For  the  case  at  hand,  these  four  vectors  associated 
d r 

with  the  corresponding  just  determined  are 


0 

' ° ' 

1 

0 

A2  = 

y 

0 

u 

0 

.a 

1 

(C.26) 


1Because  of  the  choice  of  dependent  variables 

i£ 

3P 


a2  = P / ( §£  + p * 
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and 


p[u(||)  -V] 

az  3,4 

r , dr  ^ , 3p  , dr  . 

4-vluiirp<Hr>f4 


p + [ u ( af  l 4'v]vIk 


‘ [ u ( ^7  ) ‘ v ] 

az  3,4 


8p 

aF 


(C. 27} 


where  the  last  two  vectors  satisfy  Eq.  set  (C.21)  when 

(It* 

Eqs.  (C.24)  and  ( C . 2 5)  are  substituted  for  -j±-. 

Now,  having  these  four  vectors,  the  ordinary 

dr 

differential  equations  associated  with  each  ^ are  found  by 

substituting  each  vector  into  Eq . (C.18)  and  performing  the 

inner  product.  The  differential  equation  found  for  each 

vector  is  called  the  compatibility  relation  valid  along  the 

dr  T 

associated  curve  described  by  the  particular  For  A^, 

udu  + udv  + = 0 . (C . 28) 

-V 

For  A2> 


dh  + udu  + vdv  = 0 , 


(C.29) 
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and  for  ^ , 

[U  'v;i  [udv-vdu]  - ^[u2  (g/  -w  (£)  + v2  ] 

3,4  r d23,4  dz3>4 

+ ^ [v  +u]  = 0 . (C.30) 

p az3,4 

Thus,  the  four  ordinary  differential  equations  for  the 
four  curves  defined,  by  the  ^rs  are  found. 
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LIST  OF  SYMBOLS 

Where  a symbol  is  defined  only  for  the  purpose  of 
discussion  within  the  text,  it  may  not  be  listed  here. 
Chemical  elements  are  not  listed  because  of  their  univer- 
sality. 

a Local  frozen  speed  of  sound 

0^  Specific  heat  at  constant  pressure 

C Mole  fraction  of  molecular  species 

Shape  parameter  of  mass  fraction  distribution 
E Base-bleed  parameter,  Fig.  14 

h Enthalpy 

h One-half  of  the  height  of  the  base  in  axially 

symmetric  flows;  i.e.,  (Rj  -Rj)/2;  the  base  height 
in  planar  backstep  flows 
H Total  enthalpy 

k In  the  mixing  layer,  the  mass  fraction  of  the  high- 

speed stream  in  the  local  mixture,  Eq . (16) 
k A mixing-length  constant,  Eq.  (40) 

K Equilibrium  gas  constant,  Eq . (A. 11) 

£,  A specific  x-location  in  the  mixing  layer 

I Prandtl’s  mixing  length 

m Number  of  atoms  present 

m Rate  of  mass  flow 

M Mach  number 

n Number  of  species  present 
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N 

p 

P 

P 

r 

R 

R 

Re0 

s 

t 

T 

u 

u 

V 

w 

X 

X 

y 

z 


Mole-atoms  per  unit  mass  of  mixture 
Pressure 

Recovered  pressure,  after  recompression  process 
Total  pressure 

Fictitious  partial  pressure  of  atom 
Transverse  coordinate  in  method  of  characteristics 
Radius  from  axis  of  symmetry  to  point  in  mixing 
layer 

Gas  constant  of  molecular  species 
Universal  gas  constant 

Reynolds  number  based  on  momentum  thickness 

Entropy 

Temperature 

Total  temperature 

Velocity  component  in  direction  of  developing 
mixing  layer 

Axial  component  of  velocity  in  method  of  charac- 
teristics 

Transverse  velocity  component 
Molecular  weight 

Mixing-layer  coordinate  in  the  direction  o.f 
development 

Coordinate  along  recompression  on  slipline 
Transverse  mixing-layer  coordinate 
Axial  coordinate  in  method  of  characteristics 
Number  of  moles  of  species  j [reactants) 
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3 

Y 

Y 


Number  of  moles  of  species  j (products) 

Angle  of  inviscid  boundary  with  undisturbed  free 
stream 

Isentropic  exponent 

Relative  angle  of  inviscid  boundary  with  slipline 
at  intersection  point  of  inviscid  boundaries 
Boundary-layer  thickness 
Boundary- layer  displacement  thickness 
Apparent  kinematic  viscosity 

x 

Transformed  coordinate  of  Eq . (9),  / eoF(x)/Udx 

o 

Dimensionless  mixing  variable,  ay/x 
Dimensionless  position  parameter,  o<5/x 
Boundary-layer  momentum  thickness 
Gibbs  free  energy  per  mole 
Density 

Mixing  or  spreading  parameter 
Incompressible  a,  a reference  value 
Velocity  ratio,  u/U 

Angle  between  slipline  and  axial  direction 


Subscripts 

A Axially  symmetric  value 

B Base  region 

D Dividing  streamline 

E Bleed 

i ith  species 

i Intersection  point  of  inviscid  boundaries 
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J Jet 

L Quiescent  or  base  edge  of  mixing  layer 

m Regarding  inviscid  boundary  and  coordinate 

increment  defined  in  Eq.  (15) 
s At  point  of  separation,  initiation  of  mixing 

S Stagnating  streamline 

51  Stream  1 , Fig . 1 

52  Stream  2,  Fig.  1 

U High-speed  edge  of  mixing  layer 

co  At  undisturbed  free-stream  conditions 

1 Beginning  of  recompression 

1 Of  Stream  1,  Fig.  1 

2 End  of  recompression 

2 Of  Stream  2 , Fig . 1 


no 


